Vibrational Spectrum of the B=7 Skyrme Soliton 

W. K. Baskcrville 
Department of Mathematical Sciences, University of Durham, Science Laboratories, South Road, Durham, DHl 3LE, U.K. 

(February 1, 2008) 

The finite energy vibrational normal modes of the baryon number B=7 Skyrme soliton are com- 
puted. The structure of the spectrum obtained displays considerable similarity to those previously 
calculated for baryon numbers 2, 3 and 4. All modes expected from an approximate correspondence 
between skyrmions and BPS monopoles are found to be present. However, in contrast to earlier 
calculations, they do not all have energies below the pion mass. The remaining 'breather-type' 
modes also conform to predictions, except that one predicted multiplet is not observed. 
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I. INTRODUCTION 

The normal mode spectra of multiskyrmions are be- 
lieved to play a key role in their quantisation [nl. In the 
Skyrme model Q, nuclei appear as topological solitons 
in a classical, non-linear field theory of 7r-mesons. Since 
the model is not renormalisable, a proper treatment as 
a quantum field theory is very difficult. The usual ap- 
proach to quantisation is based on the suggestion of Man- 
ton, that the low energy dynamics of n skyrmions may 
be approximated by motion on some finite dimensional 
manifold. The simplest possibility is the moduli space 
generated by the zero modes of the minimal energy so- 
lution. For B=l, collective coordinate quantisation of 
(coupled) spin and isospin gives a description of nucleons 
and the A resonance in qualitative agreement with exper- 
iment P,B. However, such a quantisation includes effects 
of order fi while ignoring effects of order %. Also, it is 
known that for real nuclei, nucleons in the ground state 
have relatively large kinetic energies. Hence, the num- 
ber of degrees of freedom quantised must at least equal 
the number of nucleon coordinates. Widely separated 
skyrmions can be labelled by their positions and relative 
isospin orientations, so a minimum of QB degrees of free- 
dom are required to describe a Skyrme configuration of 
baryon number B. However, the minimum energy always 
occurs when the skyrmions coalesce into a single large 
soliton p,p[, with a maximum of 9 zero modes: trans- 
lation, spin and isospin. Clearly, in the neighbourhood 
of such a configuration, the moduli space approximation 
must be extended to include low- lying vibrational modes. 

The first step in this direction was taken by Leese et 
al |[7|, who quantised a ten-dimensional manifold along 
the B=2 attractive channel, using an instanton approx- 
imation for the field configurations. While still not in- 
corporating quite enough degrees of freedom, their re- 
sults were significantly better than had previously been 
obtained from a simple quantisation of the eight zero 
modes. Walet [B| later extended this treatment, again 
employing the instanton approximation, to estimate the 
full normal mode spectra for B=2 and B=3. The vi- 
brational spectra for B=2, 3 and 4 have now been di- 



rectly computed by Barnes et al Wm, who were also 
able to classify the modes according to the representa- 
tions they form of the symmetry groups of the static 
solutions. A remarkable pattern emerged. The lower 
frequency modes correspond to known attractive chan- 
nel scatterings; it was also observed that these modes 
fall into representations exactly corresponding to those 
for small zero mode deformations of BPS monopoles, for 
monopole charges equal to the baryon numbers of the 
multiskyrmions. This strongly suggests an underlying 
link between the two theories, and in fact Houghton et 
al llG] have demonstrated a connection through ratio- 
nal maps. They speculate that the correspondence be- 
tween BPS monopoles and skyrmions should hold for any 
baryon number (monopole charge) , and have used this as- 
sumption to predict the lowest twenty one modes for the 
B=7 skyrmion. 

The remaining, higher frequency parts of the B=2, 3 
and 4 spectra consist of 'breathing'-type modes. The 
first mode above the pion mass is in all cases a trivial 
breather, a simple size fluctuation of the multiskyrmion. 
Above this lie more complicated breathing modes, where 
one part of the soliton expands, while another part is 
compressed. Baskerville and Michaels [Q recently pro- 
posed a simple geometric explanation of these modes, 
based on the movement of 'branch lines': lines of (ap- 
proximately) zero baryon density radiating out from the 
centre of the minimal energy solutions. This ansatz ex- 
plains all the higher frequency normal modes observed for 
B=2, 3 and 4, although it also predicts an extra triplet 
for B=4 which was not observed numerically. It is pos- 
sible that this mode, if it exists, might have a rather 
high frequency, which could explain why it was not seen. 
Baskerville and Michaels also published a complete set 
of predictions for the vibrational spectra of B=5, 6 and 
7 multiskyrmions. In general, their geometrical ansatz 
predicts 4B — 7 higher breathing modes, in addition to 
a trivial breather, for a minimal energy Skyrme solution 
of baryon number B, if B > 2. The monopole analogy of 
Houghton et al predicts a further 4B — 7 lower frequency, 
scattering-type modes. Added to the nine zero modes 
(for B > 2), this gives a total of 8B — 4 degrees of free- 



dom. This is rather more than the 6B zero modes of B 
widely separated skyrmions, throwing open the question 
as to how many degrees of freedom are really relevant in 
the vicinity of the minimal energy solution, and therefore 
how many should be quantised. It is interesting to note 
that instanton generated approximations to Skyrme field 
configurations naturally have 8B — 1 degrees of freedom. 
Houghton ||l^ recently computed the vibrational spec- 
tra of the instanton-generated 3-skyrmion, and obtained 
the same modes as predicted by a combination of the 
monopole and geometrical methods, plus an additional 
triplet, with axial vector symmetry. 

In this article, the vibrational normal modes of the 
B=7 skyrmion arc computed, using the same method as 
in ||l|,y|. Their symmetries are then classified according 
to the symmetry group of the static soliton, in this case 
the icosehedral group, the minimal energy solution hav- 
ing the shape of a perfect dodecahedron Q. The B=7 
spectrum is of particular interest since a full set of pre- 
dictions already exists for the normal modes. It thus 
provides a good test of both the monopole analogy and 
the geometrical ansatz for higher modes. In fact, all of 
the expected vibrations are observed, except for one five- 
fold degenerate multiplet. The B=7 vibrational spectrum 
may also be important for another reason. Irwin |13| re- 
cently performed a zero mode quantisation of minimal 
energy skyrmions for baryon numbers four to nine and 
seventeen. For baryon numbers four, six and eight he 
obtained ground states in agreement with experimental 
observations. However, for baryon numbers five, seven, 
nine and seventeen, the calculated ground states have 
the wrong spins. For the specific case of B=7, it may be 
possible to resolve the problem through the inclusion of 
vibrational modes. The experimentally observed ground 
state is an isodoublet with spin | , whereas the lowest al- 
lowed state in the skyrmion quantisation is an isodoublet 
with spin |. However, it is possible that a spin | rota- 
tional state could be combined with a vibrational state to 
give an allowed state. If the vibrational energy involved 
was small enough, the resultant state might have lower 
energy than the spin | isodoublet. Obviously, knowledge 
of the exact frequencies and symmetries of the low-lying 
vibrational states is a prerequisite for progress in this 
direction. 



II. COMPUTATION AND INTERPRETATION OF 
THE SPECTRUM 

The normal mode spectrum is computed using the 
same method and computer codes developed for the B=2, 
3 and 4 calculations [yj|]. Skyrme fields, being SU(2)- 
valued, lie on a 3-sphcre, and can thus be represented by 
a scalar field e SR**, with the constraint (j)°-(j)°- = 1. The 
Skyrme Lagrangian density can then be written, in terms 
of the usual dimensionless Skyrme units |17|, as 
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(1) 



where A is a Lagrange multiplier field, and cj^ — -^^ is a 
dimensionaless constant, the oscillation frequency of the 
homogeneous pion field. As in ||l|,|[, the 'standard' value 
[Q a;,r = .526 has been used. 

The Skyrme equations of motion are derived by dis- 
cretising the action S — Jd^ £ on a finite lattice, then 
varying with respect to each field component 0° (i, i, j, fc), 
at each point in the lattice. This leads in general to a 
coupled system of equations, due to the mixed time and 
space derivatives in the quartic Skyrme term. The rel- 
evant, "kinetic" part of the Lagrangian density can be 
written 
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KIN = ^ 0" K'^\d,<jy) / 
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where K""^ = 5''^ {I + {di(t)f) - di4)''d^(jJ' acts as a spa- 
tially dependent inertia tensor. It is the time dependence 
of K"'^ which leads to the coupling, with all its inherent 
numerical difficulties. However, for very small perturba- 
tions about the static solution, (j){t, x) — (?!)st(x) -I- £(t, x), 
the time dependent parts of K"^^ give rise to terms of 
order at least e'^. If e <C 1, such terms can be neglected 
to a good approximation, and the matrix K°'^ assigned 
its value at the static classical solution. The discretised 
equations of motion then reduce to the form 

r{t + l,i,3,k)^R''{t,i,j,k) + ~Xr{t,i,j,k) (3) 

where A is a constant related to the Lagrange multiplier, 
and i?" is a function of the fields at the current and 
previous timesteps only. Invoking the SU(2) constraint 
<f>1j^i(l/l_f_i — 1 yields a quadratic equation for A 



A^ + 2\R° 



R^R" -1 = 0, 



(4) 



the positive square root of which is the desired solution 
(the definitions of A and i?" can be chosen so as to ensure 
this: see [Q for details). Substituting A back into Eq. (g) 
then gives a deterministic time evolution for the Skyrme 
fields. 

The vibrational spectrum can be directly computed 
using this algorithm. The static solution is found by nu- 
merical relaxation, using fields derived from the rational 
map ansatz as a starting point |lG]. A small random 
perturbation is then applied, after which the fields are 
evolved forward for a long time at constant energy. The 
evolving fields take the form 

(j){t,x)^(j)st{x)+ Y^ e„(5„(x)cos(tj„i) + 0(e2), (5) 

modes 

where the functions i5„(x) e ^^, obeying (5„(x) • 0st(x) — 
are the normal modes, each excited with amplitude e„. 



The normal mode frequencies LUn are obtained by Fourier 
analysis of the field (/)(t,x), with respect to time, at any 
point in the box. 

Once the frequencies have been identified, the time 
evolution process is repeated. Maps of the normal modes 
(5„(x) can be constructed by performing discrete Fourier 
sums on each component of the evolving field. The space 
of perturbations has a very useful inner product 
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the matrix i^°''(x) being included in this definition to en- 
sure orthogonality between modes of different frequency: 
((5i|(52) = if wi ^ UJ2- In general, the normal modes fall 
into degenerate multiplets at discrete frequencies. Each 
such multiplet forms a matrix representation of the sym- 
metry group of the static soliton, via the inner product 
{Si\S\Sj). The symmetry operations S in general consist 
of a physical operation, for example a rotation or reflec- 
tion, which leaves the energy and baryon densities of the 
static soliton invariant, plus an isospin transformation 
of the pion fields (0^, 0^, (j)'^), which "undoes" the ef- 
fect of the physical operation on the static fields. Modes 
of the same frequency always form an irreducible rep- 
resentation (irrep) of the symmetry group, which places 
strong constraints on the allowed degeneracies. In the 
case of the B=7 multiskyrmion, the symmetry group of 
the static soliton is the icosehedral group, which has irre- 
ducible representaions of degeneracy one, three, four and 
five only. 

Determining the degeneracy and symmetry of each 
multiplet is in principle straightforward. For degener- 
ate modes, each random initial perturbation will pro- 
duce a different linear combination ^^ ei(5i(x), when one 
projects out the relevant frequency. Applying the sym- 
metries of the static soliton will yield yet other, different 
combinations. The degeneracy of a given frequency can 
be found by computing the rank of the matrix of inner 
products between different linear combinations of the de- 
generate modes. It is also possible to project any given 
linear combination out of all the others. In this way, a 
library of orthonormal modes can be built up. Once a 
full orthonormal basis has been obtained for a given fre- 
quency, it is easy to calculate its symmetry. Irreducible 
representations are uniquely labelled by the characters 
(traces J2i{^i\^\^i)) of certain symmetry operations S. 

Unfortunately, several difficulties arise in this proce- 
dure for the particular case of B=7. The main problem 
is that the simulations are carried out on a finite, cu- 
bic lattice, which breaks the dodecahedral symmetry of 
the static solution. The most immediate consequence of 
this is that some multiplets may be artificially split, for 
example a five-fold degenerate multiplet into a doublet 
and triplet of close, but not identical, frequency. The 
breaking may also cause certain modes to be generated 
more easily than others, if they fit better onto the fi- 
nite grid. Thus different initial conditions do not re- 
ally produce randomly different linear combinations of 



degenerate modes. This can effectively 'hide' some of 
the degeneracy, since the method outlined above depends 
on finding enough independent linear combinations. A 
closely related technical difficulty is that not all of the 
symmetry operations of the icosehedral group can be 
exactly realised on a cubic lattice. If the cubic lattice 
does not map onto itself under some physical operation, 
then interpolation must be used to find the new field val- 
ues. This can be done accurately enough to compute the 
character of such a symmetry, given an orthonormal ba- 
sis; but is not suitable for separating out such a basis 
in the first place. Since only 24 (out of 120) physical 
symmetries can be performed exactly on the cubic lat- 
tice, this is a severe limitation. The last (but certainly 
not least) problem is contamination between modes of 
different frequency. Forty-three vibrational modes are 
predicted for B=7, lying in at least eleven different mul- 
tiplets. These are all expected to occur within a fairly 
small range of frequencies. The finite time scale of the 
simulation leads inevitably to finite-width peaks in the 
power spectrum. Considerable overlap between different 
multiplets is therefore possible, and does in fact occur, 
as can immediately be seen from Figure (). 

The resolution of these difficulties requires a closer ex- 
amination of the icosehedral symmetry group Ih , as well 
as the particular alignment of the dodecahedral static so- 
lution on the cubic lattice. The character table for I^ is 
given in Table |[ The physical symmetries can be divided 
into ten 'conjugacy classes' of like operations. The iden- 
tity 1 and inversion i both form classes on their own. 
There are six five-fold rotational axes in a dodecahedron, 
running through the centres of opposite faces. Rotations 
by ^ about these, in either direction, form one conjugacy 
class (with 12 elements), labelled C5, while rotations by 
^ form another, C^, also with 12 elements. There are 
ten three- fold axes, through pairs of opposite corners of 
the dodecahedron. The 20 ^ rotations associated with 
these form another conjugacy class, C3. The final ro- 
tational symmetry is two-fold, about axes which bisect 
pairs of opposite edges of the dodecahedron. This gives 
15 rotations by n (there are 30 edges, since each edge 
is shared by two of the twelve pentagonal faces), form- 
ing the conjugacy class C2. The remaining symmetry 
operations are all combinations of inversion with one of 
the rotations already mentioned, and are denoted iC„, 
for rotation C„. (Note that iC2 gives pure reflections, 
in planes which contain two opposite edges, and bisect 
four pentagonal faces.) The characters for these classes 
are exactly the same as for the equivalent rotations with- 
out inversion, except that negative parity representations 
pick up a minus sign. This is because the full icosehedral 
group Ih is just the group of rotations of an icosehedron 
(dodecahedron) /, extended by inversion: /^ = / x Z2. 
The irreducible representations are labelled by their de- 
generacies and parities, where this is sufficient to dis- 
tinguish them (eg. 5~ is a five-fold degenerate multiplet 
with negative parity). The only exceptions are the four 
three-dimensional irreps, which are labelled Fi and F2, 



FIG. 1. Schematic diagram of one cube within a dodeca- 
hedron. The edges of the dodecahedron are drawn as dotted 
hnes on the faces of the cube. In reality, they would project 
slightly, and the edges of the cube would lie slightly behind 
the pentagonal faces. 



again with superscripts denoting parity. With this nota- 
tion, a vector has symmetry F^ , while an axial vector 
transforms as F^ . 

Next, we need to examine the alignment of the do- 
decahedral soliton on the cubic grid. A dodecahedron 
contains within itself five 'cubes'; that is, there are five 
sets of eight vertices, each of which form the corners of 
a cube. Each of the 20 vertices of the dodecahedron is 
a member of two such cubes. In the B=7 static soliton 
used in simulations, one of these cubes is aligned with the 
axes of the cubic lattice. See Figure (|l|) for a schematic 
diagram. It should be obvious, looking at this Figure, 
which symmetry operations map this cube to itself (and 
hence preserve the grid structure): the identity, inversion, 
three C2 rotations (about the Cartesian axes), three re- 
flections, eight C3 rotations (about the body diagonals of 
the cube), and the same eight rotations combined with 
inversion. The vital question is how to use these symme- 
tries to produce orthonormal bases for all possible mul- 



tiplets. It is very useful, at this stage, to keep in mind 
some sort of physical model for representations of various 
degeneracies. For example, the Cartesian axes will obvi- 
ously transform as a triplet (in fact F^ ) . Given a general 
mixture of all three (a vector in a random direction), the 
symmetry operations (1 + C^;), (1 + Cj,) and (1 + Cz), 
or equivalently {Cy + Cz), {Cx + Cz) and {Cx + Cy), 
will project out the x, y and z components respectively, 
where Cx is a 180° rotation about the a;-axis, etc. De- 
generacies four and five are a little trickier. The five 
cubes mentioned above form a representation, since all 
symmetries of the dodecahedron will rotate them among 
themselves. This representation is not irreducible, since 
the sum or symmetric combination of all five transforms 
trivially (as 1"*"), ie. this combination is invariant under 
all physical symmetries of the dodecahedron. However, 
the four remaining antisymmetric configurations, orthog- 
onal to the symmetric combination and to each other, do 
give an irreducible 4"*" representation. Labelling the five 
cubes {|1) • • • |5)}, a suitable orthonormal basis for the 
degeneracy 4 irrep can be written 

I^> = ^(|1> + |2)-|3)-|4)) (7a) 



|S)=2(|1> + |3)-|2)-|4)) (7b) 



|C) = -(|1) + |4)-|2)-|3)) (7c) 



1^) 



2%/5 



(4|5)-|1>-|2)-|3>-|4)). (7d) 



If the cube which is aligned with the lattice is taken to 
be 1 5), then the C2 rotations about the Cartesian axes 
act as follows: 



Cx : [A, B, C, D) — ^ (A, -S, -C, D) 



(8a) 
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TABLE I. Character table for h. t = ^%ti. 



Cy : (A, B, C, D) -^ (-A, B, -C, D) (8b) 



Cz : (A, B, C, D) -^ {~A, -B, C, D). (8c) 

Using these results, it is quite simple to construct pro- 
jection operators for the orthonormal basis 



PA^t + Cx-Cy~ Cz 



PB = 1-Cx+Cy-Cz 



Pc^l-Cx-Cy + Cz 



PD = 1 + Cx+Cy + Cz. 



(9a) 



(9b) 



(9c) 



(9d) 



It should be noted at this stage that Pa , Pb and Pc also 
act as projectors for the triplet formed by the Cartesian 
axes. What now for a five-fold basis? There are six 
C5 rotational axes, giving a 6-diniensional representation 
which decomposes in a similar way to the 5 cubes: again 
the symmetric combination transforms trivially as 1+, 
while five orthogonal antisymmetric combinations form 
b^ . LabeUing the six C5 axes {|1)---|6)}, a possible 
basis for the latter is 



|a) = -^(|l)-|2)), 



|c) = -^(|5)-|6)), 



|6) = -^(|3)-|4)), (10a) 



= -(|l) + |2)-|3)-|4)), (10b) 



2v^ 



(2|5)+2|6)-|1)-|2)-|3)-|4)). (10c) 



The action of Cx , Cy and C^ on these states is as follows 



Cx : (a, 6, c, d, e) > (a, -b, 



Uy 



,d,e) 



Cz ■ (a, 6, c, d, e) — > (—a, —6, c, d, e). 
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It is then apparent that the projection operators Pa, Pb 
and Pc defined in Eq. (||) also act as projectors for states 
I a), \b) and |c) respectively, while Pd will project out 
some combination of states \d) and |e). It may seem an 
incredible coincidence that the same symmetries project 
out the basis states for multiplets of degeneracy three, 
four and five, but in fact this is deeply related to the way 
in which the cubic lattice breaks the icosehedral sym- 
metry. Furthermore, these same projectors work for all 
degeneracy three, four and five multiplets, not just F-^^ , 
4+ and 5"*", from which they were derived. This last is 
not obvious, in fact is perhaps counter-intuitive, but nev- 
ertheless has been found to work in practice. 

The projectors defined in Eq. (m are the main tool re- 
quired to separate out the various multiplets, when con- 
tamination occurs between them. States of different par- 
ity can easily be separated, using the projectors 



P+ = t 



P_ = 1 - i. 



(12) 



However, contamination between degenerate multiplets 
of the same parity is more awkward. The coincidence 
between the projection operators means that mode \A) 
of a 4+ cannot be distinguished from mode \a) of a 5+ 
for example, or from the "x" component of F^ or F^ . 
However, it is possible to tell whether the three modes 
projected out by Pa, Pb and Pc are derived purely from 
a single multiplet, or represent some mixture, by consid- 
ering their behaviour under C3 rotations. A 120° degree 



rotation in physical space causes a cyclic permutation of 
the three projections. This is also a 120° degree rotation 
in the 'phase space' of the three modes, and since the 
modes may be aligned or anti-aligned with the Carte- 
sian axes of this space, two of the modes may pick up 
minus signs. This occurs randomly, and the chance of 
it happening the same way for two different multiplets 
which are mixed together is only ^. Hence if the overlap 
{S\Pa\S) is exactly equal to ±{S\C3Pb\S), we are almost 
certainly dealing with a pure multiplet. If, furthermore, 
{S\Pd\S) = 0, then \S) represents a pure triplet. Other- 
wise, the mode must have degeneracy four or five. The 
latter possibilities can be distinguished by the behaviour 
oi Pd\S) under C3 rotations. If the multiplet has degen- 
eracy four, then (SIPJjC'zPdIS) = 1, while for a degen- 
eracy 5 multiplet {S\PI,C3Pd\S) = —\. For degeneracy 
5, the operator (C3 — C^^)Pd ( where C3 and C^^ are 
120° rotations in opposite directions about the same axis) 
projects out a linear combination orthogonal to that pro- 
jected out by Pd, thus completing the orthonormal basis. 
The overall procedure for interpretation of the spec- 
trum is then as follows. First of all, any contamination 
by zero modes (spin or linear momentum) is projected 
out. A careful watch is kept during the simulation to 
make sure that "drift" (ie. excitation of zero modes) is 
kept within strict limits, to ensure the validity of the ba- 
sic assumption underlying the computer algorithm. Ad- 
ditionally, severe damping is applied at the edges of the 
box for the first 1000 timesteps of the simulation, be- 
fore any data collection begins. This helps to decrease 
the proportion of radiation relative to real vibrational 
modes, and also helps reduce the subsequent drift. Still 
slight contamination from zero modes occurs, which must 
be removed from the perturbations (5"(x). Next, the per- 
turbations are separated into positive and negative parity 
components. The projectors Pa to Pd are applied to the 
resulting states. The overlaps {5\Pi\5) and {5\CzPi\5) are 
calculated, as is the matrix element ((5|P^C3Pd|<^)- If 
a pure multiplet is found, an orthonormal basis for it is 
constructed. This can then be projected out of all other 
frequencies. Hopefully this will remove enough contami- 
nation that at least one other mode (of the same parity) 
will now be pure. The procedure can then be repeated, 
until there are no pure modes left. As the library of ba- 
sis states for pure multiplets is built up, a check is kept 
on how much of each frequency is accounted for. If all 
modes are suitably normalised, then 



^(<5|c.«)(c.W|<5) 



= 1 



(13) 



for all original perturbations \5), where {|a;n )} is the 
full set of orthonormal basis states. It may happen that 
some multiplets are inextricably mixed, so that it is not 
possible to extract a pure set of basis states for them. 
This will apply mainly to the "triplet" states, as modes 
\D) of 4^ and {\d), |e)} of 5* can always be separated at 



least from each other. It is still possible to tell how many 
"triplets" there are. Applying Pa, Pb and Pc to any 
one of the remaining frequencies will give basis vectors 
in the combined space of the remaining multiplets. This 
can be projected out of the other frequencies, and the 
process repeated, as for the pure modes, until all modes 
have been accounted for. In practice, this means until 
all the sums in Eq. (|l3| ) add to at least 0.995. Thus it 
is always possible to fully classify the symmetries of all 
vibrational modes in the spectrum, even if they cannot 
all be projected out cleanly. The main difference is that 
frequencies can be confidently assigned to the pure multi- 
plets, whereas mixtures will in general merely be confined 
within a certain range of frequencies. Sometimes a ten- 
tative assignment can be made if a pure \D) mode or 
{|(i), |e)} doublet appears strongly at a given frequency. 
One last technical detail should be mentioned. As was 
explained above, the symmetry operations that will be 
applied to the perturbations {|'^)}, are combinations of 
physical symmetries and isospin transformations. The 
alignment of the static solution on the cubic grid is as 
shown in Figure H. It is also necessary to determine the 
exact symmetry of the pion fields in the static solution, 
or at least, how they transform under the specific sym- 
metries needed to identify the perturbations. From the 
foregoing discussion, it should be clear that the symme- 
try operations needed to project out orthonormal bases 
for all modes are inversion, the three 7r-rotations Cx, Cy 
and Cz , plus one ^- rotation (any of the eight which can 
be realised on the cubic lattice will do) and its inverse. 
Additionally, we will need to calculate the character of 
one -?- rotation, to distinguish between the triplet rep- 
resentations Fi and ¥2- The following equations give 
the full symmetry operations used, both the physical op- 
erations on the Cartesian axes {x,y,z), and the isospin 
transformations which restore the original fields of the 
static solution. 



i : {x,y,z) — > (-x, -y, -z) 

(7ri,7r2,7r3) — > (-vri, -7r2, -tts) 
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where a = 1.0206', and 9 = 1.1410'. Note also that 
(tti, 772, TTs) are the fields (0^, 0"^, (/)^), in the notation used 
earlier in this Section. 



III. RESULTS. 

Simulations of the vibrations of the B=7 multi- 
skyrmion were carried out on a finite lattice of 65^ points, 
with side length L = 8 in Skyrme units. The resulting 
power spectrum is displayed in Figure g. 

The lowest peak in the spectrum, at frequency 0.019, 
was found to have total overlap J2i{^-Oi9\Ji){Ji\^.oi9) — 
.9911 with angular momentum zero modes. It is thus 
identified as comprising the rotational modes of the static 
soliton, which are broken from zero to finite energy by the 
discrete lattice. The same thing happened in the B=2, 3 
and 4 vibrational spectra. The next three peaks contain 
a pure 5"'" multiplet (mainly in the middle of the three 
peaks) and a pure 5^ multiplet, which is split between the 
first and third peaks, with the modes {|a), |5), |c)} mainly 
in one peak, and {|d), |e)} mainly in the other. Pictures 
of all the modes may be found in Figure |[ Both of these 
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FIG. 2. B=7 Power Spectrum. 



degeneracy 5 multiplets represent scattering type modes. 
In the positive parity 5"*" at frequency 0.302, two opposite 
pentagonal faces of the dodecahedron spht away. Pre- 
sumably these will become B=2 tori, which, if the mode 
were amplified, would move away in opposite directions, 
leaving a B=3 tetrahedron at the origin. The 5^ modes, 
split between frequencies 0.289 and 0.320, are more dif- 
ficult to decipher. It appears that one edge is pulled 
away, suggesting that the dodecahedron may be splitting 
into a single skyrmion and a B=6 soliton. Both these 
scattering processes have a "line" symmetry, which may 
explain their degeneracy. There are six natural lines in 
the dodecahedron, the C5 axes. Exciting modes in all 
six directions simultaneously however would be equiva- 
lent to the trivial breather, so that only five directions 
are really independent, hence the degeneracy. This argu- 
ment is certainly valid for the positive parity modes. It is 
not as obvious for the negative parity vibrations, but the 
similarity between the symmetries is certainly suggestive. 

The next peak, at frequency 0.358, is a pure triplet, 
with symmetry F^ . In this scattering, two pentagons 
from nearby, but not adjacent, faces of the dodecahe- 
dron are pulled away. From their relative alignments, it 
seems likely that these two tori will later coalesce into a 
single B=4 cube. This mode probably therefore repre- 
sents the breaking of the dodecahedron into two clusters, 
one with B=3 and the other B=4. A pure 4~ multiplet 
is split between the next two peaks, at frequencies 0.446 
and 0.459. This mode has a tetrahedral symmetry, with 
four corners of the dodecahedron, arranged in a tetrahe- 
dron, being pulled away. Presumably these will become 
four single skyrmions, leaving a B=3 tetrahedron in the 
centre. The opposite extreme of the motion goes through 
the dual tetrahedron, hence the negative parity. 

The next peak, at frequency 0.522, is a trivial breather 
(1+), a simple size fluctuation. After that, at 0.553, is a 
triplet of dipole breathing modes, with symmetry F^ . 
The remaining peaks all contain a mixture of several 
modes. There is a pure 4~ multiplet, spread between fre- 
quencies 0.741 to 0.842, which is a breathing mode with 
tetrahedral symmetry.The last negative parity mode is a 
pure, radiative Fj" triplet. This is the expected k — 
radiation, which because of the damping at the begin- 
ning of the run does not appear as a separate peak, but 
merely causes slight contamination (about 20%) in the 
peak at frequency 0.609. 

One more reasonably pure multiplet was identified, 
a 4+ spread across the frequency range 0.609 to 0.716. 
This is the last of the scattering modes predicted by the 
monopole analogy. It breaks the pattern previously ob- 
served in the B=2, 3 and 4 vibrational spectra, in that it 
has a higher frequency than two 'breathing' type modes. 
This may be explained by the fact that it represents a 
maximal breaking up of the B=7 soliton, splitting it into 
seven single skyrmions. This mode clearly corresponds to 
the scattering recently proposed by Singer and Sutcliffc 
pM , using the instanton approximation, as a comparison 
to the pictures in their paper will demonstrate. 



Frequency 


Degeneracy 


Symmetry 


Description 


0.019 


3 


F+ 


Broken spin zero 
modes: J+. 


0.289 - 0.320 


5 


5- 


Scattering: 1-1-6. 
One edge splits away 
as a single skyrmion. 


0.302 


5 


5+ 


Scattering: 2 + 3 + 2. 
Two opposite faces of 
the dodecahedron 
split off as B=2 tori. 


0.358 


3 


F2 


Scattering: 3 + 4. 
Two nearby tori split 
away, then recombine 
as a B=4 cube. 


0.446 - 0.459 


4 


4~ 


Scattering: 

3 + 1 + 1 + 1 + 1. 
Tetrahedral symmetry. 


0.522 


1 


1+ 


Trivial breather. 


0.553 


3 


K 


Dipole breather. 


0.609 


3 


F2 


fc = radiation. 


0.609 - 0.716 


4 


4+ 


Scattering: 1 + 1 + 1 

+ 1 + 1 + 1 + 1. 
Cubic symmetry. 


0.697 - 0.716 


5 


5+ 


Quadrupole breather. 


0.741 - 0.842 


4 


4- 


Breather. 
Tetrahedral symmetry. 


0.826 - 0.842 


4 


4+ 


Breather 

Cubic symmetry. 



TABLE II. Summary of B=7 vibrational modes. 



The remaining positive parity modes cannot be cleanly 
separated, but are a mixture of one 4+ and one 5^ mul- 
tiplet. In fact the mixing is not too severe (this can be 
seen when computing the characters), and the 5"*" modes 
can be identified as quadrupole breathing motions, while 
the 4"*" are breathers with a line symmetry similar to that 
of the 4+ scattering mode described above. 

This completes the list of all modes found below fre- 
quency 0.9. A summary is given in Table O. Definite 
frequencies are given where they can be confidently as- 
signed; where a mode is split into two or more peaks, 
or appears as part of an inseparable mixture, a range 
of frequncies is given. Modes are characterised as either 
scattering or breathing type, and an attempt is made 
to classify the former in terms of the clusters (of differ- 
ent baryon number) into which the dodecahedron breaks 
up. Some of these identifications are fairly certain, while 
others are more tentative. A full dynamical simulation 
of these attractive channels will be necessary to confirm 
the guesses made here. 



IV. CONCLUSIONS. 
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The vibrational spectrum of the B=7 dodecahedral 
Skyrme soUton has been calculated, and a full identifi- 
cation of the symmetries for all modes lying below fre- 
quency 0.9 has been made. The sums defined in Eq. jl^ ) 
were computed for all peaks, and all exceed 0.995, so 
it seems certain that nothing has been missed at these 
frequencies. 

All attractive channel scattering modes predicted by 
the monopole analogy [[lO| , were observed numerically. 
Four of the five predicted complex breathing motions 
|jll| were also observed. The 'absentee' is a predicted 
5~ breathing mode. It is possible that this multiplet 
may be lurking at a frequency higher than those consid- 
ered here. If so, it could be rather difficult to find it 
amongst the radiation which becomes more predominant 
as the frequency increases. While this mode is missing, 
the conjecture of Baskerville and Michaels ||l|], that there 
should be a total of 8B-4 modes for a multiskyrmion of 
baryon number B, cannot be said to have been conclu- 
sively proved. However, thirty eight vibrational modes 
have been identified, giving a total of forty seven modes 
on inclusion of the nine zero modes. This is definitely 
more than 6B or 6B+1 (forty two or forty three), and 
while it is possible that more vibrational modes may be 
discovered in the future, it is difficult to see how one could 
get rid of some already found. The fact that the sym- 
metries of the breathing modes observed were correctly 
predicted by ||ll[ is also persuasive. Note that the instan- 
ton inspired counting suggested by Houghton |1^] , giving 
8B-1 modes, is not ruled out either. This ansatz gives 
the same vibrational modes as the geometrical method 
of [0|, plus an extra triplet. If one multiplet may have 
been missed due to its high frequency, perhaps a second 
could be hiding up there too. The evidence would seem 
to favour the mode counting scaling roughly like 8B, at 
any rate; whether precisely as 8B-4 or 8B-1 remains to 
be decided. However, on the whole, the predictions made 
for the vibrational spectrum of the B=7 Skyrme soliton 
have proved to be accurate. 

As to the problem of the wrong spin obtained by zero 
mode quantisation of the B=7 skyrmion |13|, there cer- 
tainly exist vibrational modes of low enough frequency 
to permit hope of its resolution. In general, the more 
symmetry that is exhibited by a particular solution, the 
greater are the constraints imposed upon its possible spin 
states. Therefore, to find states with lower spin than 
the I permitted for the static solution, the vibrations 
must break its dodecahedral symmetry. The two lowest 
frequency modes found both have rather low symmetry, 
which is promising. Indeed, one can only hope that these 
vibrations still possess enough symmetry to forbid spin 
^ states. This calculation is beyond the scope of this pa- 
per, but the results presented here certainly point to it 
as a promising direction for future research. 
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FIG. 3. Contours of constant baryon density for the B=7 soliton, combined with its normal modes, as indexed by their 
frequencies in Table |l|. Where a frequency range is given in this Table, modes are labelled by the lower limit of that range. 
Note that the radiation at frequency 0.609 is not plotted, so that the mode labelled Sso9 is the 4"^ Singer-Sutcliffe scattering 
mode. All other modes are uniquely determined by their frequency or lower frequency limit, and appear in the same order as 
in Table |ll. Two different basis modes are shown for frequency 0.302, as the break-up into separate clusters shows up better in 
one, and the overall symmetry in the other. 



